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+-> Abstract 

Solutions with anti-D3 branes in a Klebanov-Strassler geometry with positive charge 
dissolved in fluxes have a certain singularity corresponding to a diverging energy density 
of the RR and NS-NS three- form fluxes. There are many hopes and arguments for and 
t-h against this singularity, and we attempt to settle the issue by examining whether this 

singularity can be cloaked by a regular event horizon. This is equivalent to the existence 
of asymptotically Klebanov-Tseytlin or Klebanov-Strassler black holes whose charge mea- 
sured at the horizon has the opposite sign to the asymptotic charge. We find that no 
such KT solution exists. Furthermore, for a large class of KS black holes we considered, 
the charge at the horizon must also have the same sign as the asymptotic charge, and 
^ is completely determined by the temperature, the number of fractional branes and the 

gaugino masses of the dual gauge theory. Our result suggests that antibrane singularities 
in backgrounds with charge in the fluxes are unphysical, which in turn raises the question 
as to whether antibranes can be used to uplift AdS vacua to deSitter ones. Our results 
also point out to a possible instability mechanism for the antibranes. 
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1 Introduction 

The backreaction of antibranes in backgrounds with charge dissolved in fluxes, like the Klebanov- 
Strassler (KS) [j] and the Klebanov-Tseytlin (KT) solutions M has been a subject of intense 
study over the past few years. A probe analysis of anti-D3 branes in the KS solution reveals 
that these branes have a nontrivial potential that drives them to polarize into NS5 branes 
wrapping a two-sphere inside the large three-sphere at the KS tip [3] . On the other hand, if one 
tries to go beyond the probe approximation and obtain the backreacted solution corresponding 
to smeared anti-D3 branes at tip of the KS solution, a surprise awaits: both the first-order 
backreacted solution [2jj7], as well as the fully backreacted solution [8] have a singularity. 

The fate of this singularity is crucial. Adding anti-D3 branes to a KS-like throat is the most 
generic way to uplift the vacuum energy of the AdS vacua that come of out of string theory 
flux compactifications with stabilized moduli [9] and obtain deSitter vacua. If the singularity of 
the anti-D3 solution is not physical, this implies that anti-D3 branes cannot be used to give KS 
metastable vacua and to uplift the AdS vacua to dS. Since the other known uplift mechanisms 



(such as F or D-term uplifting 10,11 or Kahler uplifting 12,13]) are much less generic, this 
would imply that string theory does not have a landscape of dS vacua. Thus the fate of this 
singularity is closely intertwined with that of the landscape. 

In both the first-order and the fully backreacted solution, this singularity comes from three- 
form RR and NS-NS field strengths whose energy densities diverge. There have been quite 
a few arguments both in favor and against this singularity. The arguments in favor of this 
singularity |14j are based on the self-consistency of the probe approximation of j3] and on the 
fact that the divergent energy of the singularity has a finite integraQ The arguments against 
it are that the self-consistence of the probe approximation does not imply the existence of 



metastable vacua when backreaction is taken into account 15,16]. Furthermore, the finiteness 
of the integral of the divergent energy density near a singularity can hardly constitute a criterion 
for accepting it: the negative-mass Schwarzschild solution also has a singularity with finite 
energy, and yet has to be discarded as unphysical | IT PI 



1 In 14 it was also argued that the singularity may be an artifact of first-order backreaction, but this has 



been shown not to happen |8j. 

furthermore, there are similar singularities near anti-M2 18 19 and anti-D2 branes 20 , and for those 
singularities both the energy density and its integral diverge. 
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It has also been argued that this singularity signals the tendency of the branes to polarize 
(as it happens in the probe approximation |3j), and one can therefore hope that this singularity 
could be resolved by brane polarization |21| a la Polchinski-Strassler 22 . However, the recent 
calculation of |23|^] pours cold water on this hope: neither the smeared anti-D3 branes, nor 
the localized ones polarize into D5 branes wrapping the S 2 of the warped deformed conifold 
at a finite distance away from the KS tip despite the fact that all the terms of the Polchinski- 
Strassler polarization potential are there. Since in Polchinski-Strassler there are always multiple 
channels for brane polarization, the absence of a D5 brane polarization channel for localized 
branes suggests that the NS5 channel found in a brane probe approximation in [3] is also not 
present in the backreacted solution. 

Given that all the calculations made so far that could have either confirmed or invalidated 
the arguments in favor of this singularity have given negative results, the main hope that is left 
is that some hither unknown physical phenomenon will come to its rescue and resolve it. Even 
if this "resolution by mystery" proposal has not yet been articulated, it does appear to us that 
disproving it beforehand might once and for all settle the discussion. 

The clearest argument that a given unknown singularity can be physical has been formulated 
by Gubser in 28 , who argued that if the singularity can be cloaked by an event horizon it is 



physical, and conversely if a singularity cannot be obtained by "turning off" a black hole horizon 
than it is not physical. Hence, if the anti-D3 singularity were physical, one would expect that 
there should exist a black hole in Klebanov-Strassler and/or Klebanov-Tseytlin whose charge 
has the opposite sign to the charge dissolved in the fluxes. 

Our purpose is to revisit the KT and KS black holes constructed numerically in the past and 
to show that no such black hole exists. The existence of Klebanov-Tseytlin black hole solution 
was first proposed in 30 and then was constructed in perturbation theory around the black hole 
in Klebanov-Witten 29 in 31 32 . However, the equations underlying the nonlinear solution 
cannot be solved analytically, and the full solution was constructed numerically in 35 37 
Interestingly-enough, extending the ansatz to describe black holes in Klebanov-Strassler is not 
so hard, but it turns out that these solutions (which would be dual to a deconfined phase with 
spontaneously-broken chiral symmetry) do not exist [40]! The only way to build a black hole in 
Klebanov-Strassler is to turn on one or two non-normalizable modes corresponding to gaugino 
masses in the dual theory, which break explicitly the chiral symmetry 
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In this paper we review the one-parameter family of KT black holes, and the three-parameter 
family of mass-deformed KS black holes and we calculate the D3 charge at the horizon. We find 
that this charge does not have an opposite sign to the asymptotic charge of the solution. In fact, 
for a KT black hole of a given temperature the value of this charge is not a free parameter, as 
one might expect naively from the perturbative analysis of 32 . Hence, if one imagines keeping 



the temperature fixed and lowering probe charges into the black hole, the configuration will 
settle back to the original charge. This comes essentially because the KT solution has charge 
dissolved in the fluxes, and an over-charged or an under-charged black hole can expel or absorb 
charge from the surrounding charge in the fluxes to bring back its charge to its initial value. 
A similar story happens for the mass-deformed KS black hole: if one adds positive or negative 
charge to the black hole keeping the temperature and gaugino masses fixed, the black hole will 
interact with the surrounding flux and return to its original charge. 

In section ^ we discuss the KT black hole and examine the symmetries of the background, 
the single-parameter nature of the solution and the relationship between its temperature and 
the charge. In section (|3| we explore in more detail the three-parameter mass-deformed KS 



24 



27 



3 As well as the earlier analysis of anti-D6 singularities 
4 Other attempts to construct such black holes were discuss"ecT in 33 34 . 
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black hole of 40 and explain the relation between gaugino masses, temperature and charge. 



We conclude with a few comments and suggestions for future work in section Q. 

2 The Klebanov-Tseytlin black hole 

2.1 The cascading gauge theory and its chiral symmetric phases 

The M — 1 supersymmetric SU (N) x SU(N + M) cascading gauge theory can be realized 
as the world-volume theory on N regular D3-branes to which we add M fractional D3-branes 
(wrapped D5-branes) to the apex of the conifold singularity. At low temperatures this cascading 
gauge theory spontaneously breaks a discrete chiral symmetry, which corresponds in the dual 
bulk to the deformation of the conifold singularity. The full warped deformed conifold solution 
dual to this cascading gauge theory was constructed by Klebanov and Strassler (KS) in [l]. 
One can also construct a singular solution dual to the chirally-symmetric phase of this theory 
(the Klebanov-Tseytlin (KT) solution |2|). The two solutions have D3 brane charges dissolved 
in the fluxes. 

This cascading gauge theory is an ideal testing ground for understanding the deconfinement 
phase transition of strongly-coupled QCD-like gauge theories. The KS solution is holographi- 
cally dual to the confined chiral-symmetry-broken phase, and the high temperature deconfined 



phase with unbroken chiral symmetry was argued in 30 to be dual to a black hole added to 
the KT solution. This black hole was constructed in a perturbative expansion [3lJ[32j[35j[37] in 
the fractional brane number M around the black hole in Klebanov- Witten |29l and then at full 



nonlinear level in 36 37 using numerical methods. 

We will review this solution with a highlight given to its properties that are relevant for 
our study: we emphasize that this is a one-parameter of solutions and we show that this black 
hole will always have positive Maxwell D3-brane charge at the horizon. Hence, there are no 
KT black hole solutions with anti-D3 brane charge, contrary to what one might have expected 
from the perturbative story. 

2.2 The Klebanov-Tseytlin black hole 

We review the construction of the Klebanov-Tseytlin black hole, which is a solution of 10- 
dimensional type IIB supergravity that asymptotes to the Klebanov-Tseytlin (KT) background 
[2], and is holographically dual to a high-temperature deconfined chirally-symmetric phase of 
the cascading gauge theory. Recall that the IIB supergravity field content includes, besides the 
gravitational field, a dilaton $ = Ing, a NS-NS flux H3 = dB 2 , and R-R fluxes F3 and F5 (the 
axion vanishes, Cq = 0). These solutions have a Z 2 m chiral symmetry, a U(1)b symmetry and 
a SU (2) x SU (2) global symmetry. In the absence of the black hole, the KT background is a 
direct product of .M5 with metric and a squashed T ' with the radii of curvature of 
and the fluxes if 3, F5 varying logarithmically in the radial coordinate. 

The most general ansatz (in the Einstein-frame) that is a deformation of the KT background 
(2) which preserves its SU(2) x SU(2) x Z 2 p x U(1) B symmetry, is 35 36 



ds 2 w = g^(y)dy»dy u + ds 2 Thl , ds 2 TlA = Q\{y)e% + Q 2 2 (y) J] (e 2 a + eJJ , (2.1) 
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and (we set a' = 1) 

K(y) 

F 3 = P e^A (e 01 A e 01 - eg 2 A e 02 ) , B 2 = -^p~ (e* A e 01 - e e2 A e H ) , 2 ^ 
-P5 = ^5 + ^5 , ^5 = ~K(y) A e 01 A A e 02 A , 

where y denotes the coordinates of Ai 5 (greek indices /i, z/ run from to 4) and ds^i,! is the 
line element of the warp-squashed T 1,1 , with the one-forms e$, eg a , e$ a [a = 1, 2) given by: 



e ^ = \ f# + ^ cos fl d0 a j , e<j a = -^d9 a , e 0a = sin6» a d0 a . (2.3) 



The range of the T 1,1 coordinates is < ip < Att, < 9 a < it and < <p a < 2tt. The dimensionful 



constant P in (2.2) is related to the quantized dimensionless units of flux M entering in the 



rank of the gauge groups of the dual field theory. We work in the normalization where 

where G5 = Gio/vol r i,i is the 5-dimensional Newton's constant obtained after the dimensional 
reduction on T 1,1 done next. 

With this ansatz, we can do a standard Kaluza-Klein reduction of the type IIB action to 



five dimensions and get the effective action 35 



1 K 



(2-5) 



2 ^^^2 

where Rio is the 10-dimensional Ricci scalar related to the 5-dimensional Ricci scalar R$ by 



R w = R 5 - 2nj 1 g Xl/ [ VaV^i - 8fiJ Y"( v A v,n 2 



(2.6) 



2 ,1 o2 o-4 
2 ) 



+ 24 fi 2 2 - 4 fif Q 
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The associated equations of motion are 
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-g 



2P 2 



^gg» v d v K 



K 








2^* 



n ie -*(dK) 2 P 2 e 



:-g ' L J 4P 2 

n^Rs - 12Ql(dn 2 ) 2 + 24fi 2 - 12Q 2 - 8tt 



i p 2 p* 



(<9if) S 



+ 



if 2 



4P 2 2fi 2 ^ 



4fii^i? 5 - 8^Qfii - 24fiifi 2 Qft 2 - 24fi 2 «9fii<9fi 2 - 24fiifi 2 (dfi 2 ) 2 + 48fiifi 2 

5 5 

2/T 2 

-2Q 1 ^(9$) 2 1 



o — ^ 1 



(2.7) 



where {OF) 2 denotes g^ u d^Fd u F and D 5 is the Laplacian in .M5. 

An ansatz for the A4.$ gravitational field that is tailored to search for the KT black hole 
with horizon located at x — 1 and asymptotic KT background at x — > is |36] 



r/<? 2 



(2.8) 



/i" 1/2 (2x - x 2 )- 1 ' 2 [-(1 - a;) 2 cft 2 + dx\ + dx^ + dxl] + G^dx 2 , 

and it is also useful to rewrite the squashed T 1 ' 1 factors as 

C^/^V^, ^2 = h l ' A ^h. (2.9) 

In these expressions, ft, / 2 , ^3, G xa . are functions of the compact radial coordinate x. Recall 
that the solution is fully specified once we find in addition the dilaton $ = lng(x), and flux 



forms (2.2), which are determined by a single function K(x). 



In these conditions, it follows from the equations of motion (2.7) that G xx is given by an 



algebraic relation that is a function of {h, / 2 , ^3, K, g} and their first derivatives. So, in total, 
to determine our black hole we need to solve a system of five second-order coupled and non- 
linear ODEs for {h, / 2 , / 3 , K, g}. The explicit algebraic relation for G xx and the system of 



five equations of motion are explicitly written in equations (2.6)-(2.12) of 36 . An important 



observation is that each of these five equations is second order. Therefore the total differential 
order of the system is 10. 

To find the KT black hole we solve the equations of motion subject to appropriate boundary 
conditions at the horizon, x — > 1, and at the asymptotic boundary x — > 0. 

The IR boundary condition (horizon) is quite standard in black hole physics. Formally, 
we rewrite our ansatz in ingoing Eddington-Finkelstein coordinates (appropriate to extend the 
analysis through the horizon) and we require regularity of all our fields at the horizon in this 
coordinate system. This amounts to say that all our functions must have a series expansion in 
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(1-x) 2 (36), 



x 



2n 



n=0 

oo 

^=E^(i-^) 2n 



\2n 



.r 



, 2/i 



?1=0 



n=0 

oo 

n=0 



n=0 



(2.10) 



\2n 



with the boundary condition simply requiring that the coefficients {Hq, aft, bft, k^, g^} are con- 
stants. Solving the five equations of motion perturbatively around the horizon up to arbitrary 



order one finds that only a few of the coefficients in (2.10) are independent, with all the others 



being a function of these. Concretely, one chooses the 6 IR independent parameters to be 
IR independent parameters (6): {hft, aft, bft, A;q , 

9o, «?}■ (2-11) 

Consider now the UV asymptotic structure at i — > 0. The UV boundary condition is 
straightforward: we want the KT black hole to approach asymptotically the KT solution j2]. 
The latter is the leading order contribution of a power series expansion in x and ln(x) of the 
equations of motion 



36 



h —hn — 



P 2 9o 



oo n—1 



oo n—1 



in(x) + E x " /2 lnk ^> /2 = a o+EE a ^ %n ' 2 lnfc ( x )' 



71=1 k=l 



oo n—1 



oo n—1 



n=l k=l n=l k=l 

1 oo n—1 

K =Ah Q al - -P 2 g - -P 2 g ln(x) + ^ ^ K n>k x n ' 2 ln fe (x). 



n=l fc=l 
k, 



n=l fc=l 



(2.12) 



The presence of the ln x terms in h, K in the KT leading contribution is responsible for the 
logarithmic running of the fluxes and M.5 radii of curvature. 



As expected, not all coefficients in the expansion (2.12) are independent. There are 4 in- 



dependent asymptotic parameters already present in the KT solution and that we choose to 



be {P, g , a , h }: P is related to the quantized flux M by (2.4), g Q is related to the dimen- 
sionless parameter of the cascading gauge theory, and there are two combinations of the other 
2 parameters that are related to the temperature and dynamical scale of the cascading gauge 
theory. More concretely, one has = 4:7rG^ sT, where s and T are the entropy density and 
temperature of the solution, which can be computed directly from the horizon area and surface 



gravity using the expansion (2.10). In addition we can replace the independent parameter ho 



by a new dimensionless parameter k s defined as, 



4fr ao 
P 2 9o 



1 



In 



AttG 5 sT 
A 4 



(2.13) 



This relation defines the dynamical scale A of the cascading theory. The other UV independent 
coefficients correspond to vacuum expectation values (VEVs) of the operators dual to the fields 
we are solving for. After using the conformal anomaly equation [35], one finds that there are 
four VEV independent parameters, chosen to be {02,0, 92,0, a 3,o, 04,0}- 

We conclude that there are eight UV independent parameters, that determine the coefficients 
in (2.12) to any order we wish. These can be chosen to be {P, g , a , k s , a 2 ,o> 92,0, a 3,(b 0-4,0}- 
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We can however use the symmetries of the problem to get rid of some of these parameters. 
First observe that that the ansatz (2.8) is invariant under the scaling symmetry^] 



(t, x, x) — > (A 2 t,x,\ 2 x) 



h^\~ 2 h. 



/2,3 — > A/2,3 



g^g K K . (2.14) 



We use this scaling symmetry to set ao = 1. In addition, we are working in the supergravity 
approximation (where we neglect g s and a' corrections). In this approximation, the action and 
the associated equations of motion do not depend separately on P 2 and g but only on the 
combination P 2 g^. We can thus set go = 1. Furthermore, when we neglect a' corrections, the 
action is multiplied by a constant when we rescale the ten dimensional metric by a constant 
factor (and rescale the p-forms accordingly), so that the equations of motion are left invariant. 
This transformation acts on our variables as 



h ->■ \- 2 h. 



/2,3 — > A 2 / 2 ,: 



K -»■ A 2 K 



9^9 



P^XP 



(2.15) 



As long as we work in the supergravity approximation we can thus use this symmetry to set 



P = 1 and use the scaling (2.15) to generate the solutions for any other value of P. 



To sum, after using the symmetries of our system we find that we are left with 5 independent 
UV parameters, 



UV independent parameters (5): { k s , 02,0, #2,0, 03,05 a 4,o}- 



(2.16) 



At this point we can ask how many parameters we need to describe the KT black hole 
solution. We have a total of 6 + 5 = 11 IR and UV independent parameters (2.11) and (2.16), 



and the equations of motion are a system of total differential order 10. We thus conclude that 
the KT black hole is a one-parameter family of solutions. In the numerical search of the KT 
black hole we will take this dimensionless parameter to be k s . Equation (2.13) can then be 



used to translate the results we obtain in terms of the dimensionless temperature T/A. 

The black hole is constructed numerically using a standard shooting method. One fixes 
the microscopic free parameter of the theory, k s , and wants to find the other ten IR and UV 
independent parameters (2.11) and (2.16). The equations of motion have two critical points: 
at the horizon, x — 1, and at the boundary, x — 0. Consider first the horizon. Using the IR 



series expansion (2.10), one constructs the solution in the near-horizon region up to the tenth 
order (n = 5 inclusive) in the radial distance to the horizon. This solution depends on the 6 



IR parameters (2.11 ). One then integrates numerically the six radial second order ODEs, using 



a Runge-Kutta method, up to a large radial distance. The procedure is repeated, this time 
at the asymptotic critical point where one starts by obtaining the asymptotic solution up to 
ninth order (n = 9, inclusive) using the UV series expansion (2.12). Fixing k s this asymptotic 
solution is a function of the four UV parameters {a 2t o, 92,0, 03,0 , a 4,o}- One integrates this 
solution down to very small values of the radial distance. In the overlapping region of the two 
solutions we then do their matching, more precisely at x — 0.5. The requirement that both 
the set of five functions {h, f'2, fs, K, g} and their first derivatives must be continuous gives 
ten conditions that we use to fix the 10 IR and UV independent parameters (for a given k s ). 
The whole process is now repeated for other values of k s to generate the results that will be 
presented in the next subsection. 



5 This scaling symmetry leaves ZioOo invariant. This means that ft-o^o is a function of the dimensionless 
parameter of our theory, which is the ratio between the temperature and the dynamical scale A. This motivates 



the definition (2.13) 
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2.3 Properties of the Klebanov-Tseytlin black hole 

We are ready to discuss some physical properties of the Klebanov-Tseytlin black hole that are 



relevant for our study. For a detailed account of other properties see 36 37 . 

We present the results as a function of either the microscopic parameter k s (which defines 
uniquely the solution) or as a function of the dimensionless quantity T/A, where T is the 



temperature of the solution and A is the dynamical scale of the theory, related to k s via (2.13). 
As explained above, the numerical construction of the solution exploits the scaling symmetries 
of the theory, and is done at P = go = a = 1. One can then restore the correct powers of P 



using the scaling symmetry (2.15), which also introduces a factor of go/a together with every 



factor of P 2 . In order to relax the a$ = 1 condition we then use (2.13) which implies that all 
the dimensionful quantities must be computed in units of A = e~ hs ^ 2 . 

From the surface gravity we can compute the black hole temperature, which in A units is 

T e^ 2 




(2.17) 



and enables an immediate translation of a dependence on k s into a dependence on the tem- 
perature. The entropy density is computed from the horizon area. The energy density £ and 
pressure V of the KT deconfined plasma phase can be read from the expectation values of 
the appropriate components of the holographic stress-tensor T ab : £ = (T tt ) and V = (T x . x .). 
Moreover, since there are no chemical potentials, the free energy density is a simply J 7 = —V. 
These quantities are dimensionless if expressed in A units, and given by 



47rG 5 s fT\ 4 2ks 32tt 4 £ 3/ 4a 2 , \ 2ks 32tt 4 J 7 3 fa 2>0 7 



P 2 g 2 T 3 \AJ " ' 81 A 4 4V 1 + 7a ; eS ' 81 A 4 7 \ a 12 1 ' " 

(2-18) 

We can now describe the KT black hole using these thermodynamic quantities. The mi- 
croscopic parameter k s has no positive upper bound. For large positive values of k s , the tem- 
perature is large and then decreases as k s decreases (see Left panel of Fig. [IJ. However, as we 
keep decreasing k s , the temperature reaches a minimum at k s = fc" = —1.230(3) (green dot in 
Fig. [T]) and, for even smaller k s , it increases until a critical lower bound, k s = kf %t , is reached 
(represented by the last numerical red point in Fig. [TJ. As this critical value is approached the 
curvature invariant (like the Kretschmann scalar in the Einstein frame) grows without bound, 
as shown in the Right panel of Fig. [T], which implies that the solution will develop a naked 
singularity. Below k s = kf lt there are no longer KT black holes. In Fig. [TJ the inset plot shows 
the behavior of the temperature around its minimum value. In all the plots of this Section, the 
numerical red dots describe the region with k s < k u and the blue dots correspond to k s > k u . 

The Left panel of Fig. [2] shows what happens to the energy density as a function of the 
temperature in the most interesting region of the parameter space. For larger temperatures it 
grows monotonically. On the other hand, the free energy density as a function of the tempera- 
ture is plotted in the Right panel of Fig. [2] For large values of the temperature the KT black 
hole solution is the preferred thermodynamic phase since its free energy is negative. However, 
it vanishes at 

T c = 0.6141111(3) A, (2.19) 



6 To pinpoint with higher accuracy the critical temperatures of the green (T = T u ) and magenta points 
(T = T c ) we have made runs in the code in a small window around these points with much higher resolution 
than the one presented in the figures. 
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35 000 




Figure 1: Left Panel: The dimensionless temperature ? as a function of the microscopic 
parameter k s . The inset plot zooms it in the neighbourhood of T = T u (green point). The red 
dots correspond to k s < k u and the blue dots correspond to k s > k u . To the left of the last 
red dot there are no black holes. Right Panel: The Kretschmann scalar at the horizon as a 
function of the microscopic parameter k s . 



and is positive for smaller temperatures: a first order deconfinement / confinement phase tran- 
sition occurs at T = T c (represented by the magenta point in Fig. [2J. The inset plot in this 
right panel shows how the free energy behaves in the vicinity of T = T u (green point). We see 
that in this region the free energy keeps positive but precisely at T = T u (kg) with 

T u = 0.8749(0) T c (2.20) 

we find a cusp, and a continuous phase transition occurs. To the right of this cusp, and 
within the temperature window T u < T < T(kf" lt ), we have two phases. The red phase, with 
fccnt has always higher free energy for a given temperature. So it is never the 

thermodynamically preferred phase in the canonical ensemble. Note also that the solution is 
regular at the cusp (in particular, in the equivalent plot of free energy versus k s , this cusp 
corresponds to a maximum inflection point of the curve). 

Plotting the dimensionless entropy density ^ as a function of the dimensionless energy 

density ^ff-xir we find that the two branches of solution can never have the same energy and 
entropy densities (see Left Panel of Fig. |2|. 

So, why is T = T u so special? The answer relies on the study of speed of sound waves in the 



cascading plasma 37 . In the gravitational description, this speed of sound c s can be computed 
from the lowest quasinormal mode dispersion relation in the sound channel. Alternatively, in 
the dual gauge theory description, it is computed from the standard thermodynamic relation 

It turns out to be 

1 7_i2^M_6P 2 ^p 

c 2 _ _ ao dk a ?2 21) 

s ~ 3 7 + 4 I da 2,o ' v ■ ) 

ao dk s 

In the vicinity of T u (k%) (determined by the condition c s = 0) we find that c s ~ (k s — and 



c ~ 



(l — ^f) 1,/4 . Thus, the speed of sound becomes imaginary for k s < k^ (red dots), as shown 
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Figure 2: Left Panel: The dimensionless energy density as a function of the dimensionless 
temperature. The inset plot shows this energy in the vicinity of the onset of the perturbative 
instability of the cascading plasma (where T = T u , the green point). Right Panel: The 
dimensionless free energy density as a function of the temperature for temperatures at the 
deconfmement transition (T = T c ; magenta point) and around it. The inset plot shows the free 
energy density in the vicinity of the perturbative instability of the cascading plasma (T = T u ; 
green point) and all the way down to the last red point that represents the last KT black hole 
that exists. In these figures, the red/blue dots have the same meaning as in Fig. [I] 



in the right panel of Figure |3| Consequently, the specific heat cy = ^ of the cascading plasma 
diverges at the threshold of trie linear instability, T = T u . We conclude that the KT deconfined 
plasma phase becomes unstable at this temperature. Note that from the inset plots of Fig. [2] 
one has that §^| T T is finite while §§| T T diverges. This divergency in the change of the 
energy density is responsible for the vanishing of the speed of sound. 

To summarize, the KT black hole is dual to the chirally-symmetric deconfined phase of the 
cascading gauge theory. This is the thermodynamically dominant phase of the theory at high 
temperatures. As its temperature decreases, we reach a critical point where its pressure (free 



energy) first vanishes and then becomes negative (positive). At this critical temperature (2.19), 
a first order deconfinement /confinement phase transition occurs. This is a non-perturbative 
phase transition since it proceeds via the nucleation of bubbles of the confined phase. Below 
this temperature, the KT deconfined chiral symmetric phase is still a metastable phase of the 
system all the way down to T u = 0.8749(0)T C , where it joins a perturbatively unstable branch 
(red dots in the plots) of the theory with negative specific heat. 

There is another property of the black hole that is particularly important for our study. 
The KT black hole is also characterized by its D5-brane and D3-brane charges. The (Maxwell) 
D5-brane charge, also called the fractional D3-brane charge, is simply 

2 3 / 4 7T 

Qm x = M= — m P, (2.22) 
3G 5 



where we used (2.4). The definition of the D3-brane charge is more subtle since the logarithm 



dependence of the fields makes a naive definition diverge, and we need a suitable regularization 
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Figure 3: Left Panel: The dimensionless entropy density as a function of the dimensionless 
energy density in the vicinity of T = T u . Right Panel: Square of speed of sound c 2 s in the 
vicinity of T = T u . 



procedure. For a supergravity solution with non-trivial Wess-Zumino terms, one can define 
three different types of charges [38j[39], that we now discuss. For a solution asymptoting to the 
KT background and with a horizon, the dimensionless D3-Page charge is [a 1 = 1) [6] 

where the integration is done over the T 1,1 . The D3-Page charge is conserved and is independent 
of the radius at which it is evaluated. 

We can also define the Maxwell D3-charge of the solution as 

<&- = (dW^*- (2 ' 24) 

where the integral is again performed over T 1,1 , but this time with the integrand evaluated at 
a certain cut-off x = x c . There are two physically distinct contributions to the Maxwell charge, 
from the black hole (or from mobile branes when no horizon is present), Qf,, and from charge 
dissolved in the fluxes (Qf): 



(2.25) 

, A 2 \2 I H 3 AF 3 . 

(4 7T ) J Me 

where means that the integrand is evaluated at the horizon x = 1, and Ai^ is the "bulk" 

spacetime spanned by the coordinates on T 1 ' 1 plus the radial coordinate. The Maxwell charge 
depends on the scale at which it is measured, but if we fix a holographic screen, we expect 



QMax 
D3 

Qb 3 



3 + Qf , 



(4tt 



2^2 



l.i 



Qf 



D3 



(4tt 



2\2 
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physical processes to preserve its value at the screen. In particular, for a given scale, it must 
be the same if two solutions are to describe different vacua of the same theory. 

For the KT black hole, we compute the Maxwell charge at the black hole horizon using 



(2.10), to obtain: 



Q 



D3 



27n 



k'n 



(2.26) 



As described previously, the value of fcg is determined numerically and it is always positive, for 
all values of the single parameter k s that parametrizes the KT black hole family. Therefore, 
the KT black hole can only have positive Maxwell D3-brane charge at the horizon, as explicitly 
shown in Fig. [4j In particular, the D3 charge at the horizon Qf 3 can never have opposite sign 
to the asymptotic D3 charge Q^ x ■ 
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Figure 4: The dimensionless D3-brane charge of the black hole Qf 3 as a function of the 
dimensionless temperature. 



3 The mass-deformed Klebanov-Strassler black hole 

3.1 Broken-chiral-symmetry phases of the theory. Mass-deformed 
cascading theories. 

So far we have discussed just deconfined/confined phases of the cascading gauge theory with 
unbroken chiral symmetry. From the holographic gauge theory perspective, the chiral symmetry 
preservation is explicitly manifest in the absence of expectation values for dimension-3 operators 
in the dual thermal states. On the other hand, one knows that at zero temperature the cascading 
gauge theory confines into a chiral-symmetry-breaking (xSB) phase, with the supergravity dual 
of this phase being the Klebanov-Strassler (KS) warped deformed conifold solution [l]. In the 
UV region, this solution asymptotes to the KT solution but, deep into the IR bulk, the Z 2 m 
chiral symmetry is spontaneously broken to Z 2 (by gaugino condensation in the pure SYM 
limit of the theory) and the KS solution also breaks the U(1)b symmetry, while preserving 
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the SU(2) x SU(2) symmetry. If we heat the theory, the confined phase persists: we have 
now a thermal gas of hadrons described by the KS solution with a thermal identification. A 
question that immediately emerges is then whether this theory also allows for the existence of 
a finite-temperature deconfined phase, this time with broken chiral symmetry, that we would 
naturally call a KS black hole. This question does not come alone. A priori there is no reason 
for this confinement /deconfinement phase transition of the broken chiral symmetry phases to 
occur at the same critical temperature T c — see (2.19) — of the phase transition for the chiral 
symmetric phases. 

In supersymmetric gauge theories chiral symmetry breaking necessarily occurs at temper- 
ature at least as high as the confinement temperature. The question remains concerning the 
thermal competition of the two deconfined phases: what is the critical temperature, T x sb, at 
which the KT black hole would exchanges dominance in the partition function with the KS 
black hole? That is, what would be the relation between the critical temperatures for the 
confinement/deconfinement and the broken- /unbroken-chiral-symmetry phase transitions? 

These questions were analyzed in detail in |40], whose conclusions we summarize next. We 
will soon realize that some of the above-described intuition needs to be re-evaluated. Consider 
first the KT deconfined phase with unbroken chiral symmetry. A good way to learn whether 
this phase can condensate into a broken chiral symmetric phase is to study chiral- symmetry- 
breaking (xSB) deformations of this phase. In the supergravity description this amounts to 
study linearized gravitational perturbations about the KT black hole solution that break the 
chiral symmetry of the background. To make our life easier we want to restrict to the simplest 
sector of perturbations that break the desired symmetries while keeping frozen the fluctuations 
of fields irrelevant for the discussion. Ref. [40] found this sector of deformations that solves the 
linearized equations of motion of the system. The boundary conditions further constrain the 
quasinormal mode spectrum of frequencies. The key observation is that the imaginary part of 
this frequency spectrum changes sign at a critical temperature 

T xS B = 0.882503(0)T c >r u . (3.1) 

Actually this x^B instability is a Gregory-Laflamme type of instability since it requires the 
breaking of the translational invariance along the planar spatial directions transverse to T 1 ' 1 . 
[] We seem to have all the ingredients to expect a new branch of black holes in a phase 
diagram of solutions of the theory. That is, at the threshold of x^B instability, we expect a 
natural merger of the KT black hole with what would be the KS black hole. This is where 



the above intuition proves to be dramatically wrong. Indeed, Ref. 40 numerically searched for 
gravitational solutions describing the homogeneous and isotropic states of the cascading plasma 
with spontaneously broken chiral symmetry. The idea was to construct the KS black hole by 
deforming the chirally symmetric KT black hole for T < T x sb along the tachyonic directions 
revealed in the linearized analysis. This attempt was unsuccessful: there are no KS black hole 
solutions in the cascading gauge theory. 

At this point, we still need to add a twist to the story. So far, we have been silently 
assuming that the KS deconfined phase we were searching for was breaking the chiral symmetry 
spontaneously. If we are less demanding, we can still search for thermal deconfined phases of 
the mass-deformed cascading gauge theory that explicitly break the chiral symmetry. To get 
the mass-deformed cascading gauge theory one introduces the mass terms 

Tfl ' 



7 As an interesting side note, recall that the deconfined KT plasma is thermodynamically stable down to T u 
and the fact that the x^B Gregory-Laflamme instability occurs at T x sb > T u provides an explicit example of a 
violation of the Gubser-Mitra conjecture. 
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for the gauginos (the J\f = 1 fermionic superpartners of SU(N + M) x SU(N) gauge bosons). 
These mass terms explicitly break both the supersymmetry and the chiral symmetry, and the 



theory does have a homogeneous and isotropic deconfined thermal phase. The mass terms (3.2) 



are the couplings of the two dimension-3 operators that explicitly break the chiral symmetry 

Oi = Oi(^), j = 1,2, (3.3) 
and, in the chiral limit /ij — > 0, the expectation values for the condensates vanish as well: 

lim 0foi,) = O. (3.4) 

In the supergravity description, the homogeneous and isotropic deconfined broken chiral phase 
is described by what we call the mass-deformed KS black hole. This black holes was numerically 



constructed at the full non-linear level in 40 . It is reassuring that the same numerical code 
that does not find the KS black hole dual to what would be a spontaneously broken chiral 
symmetric deconfined phase of the theory (/i; = 0), does find the mass-deformed KS black hole 
(/ij 7^ 0). This definitely guarantees that the statement that KS black holes do not exist is not 
a consequence of a somehow incomplete numerical search. 

We can now turn back to the linearized study of the xSB physical excitations of the KT black 
hole. The original study of |40] included linearized perturbations with the mass-deformations 
turned-on (these mass-deformations are dimension-3 operators and can be read from the coeffi- 
cients of the non-normalizable modes that decay asymptotically as 1/r; their normalizable part- 



ners give the VEVs for these dimension-3 operators (3.3)). This study finds that at T = T x sb 
the KT black hole becomes linearly unstable. As described below, the threshold of this insta- 
bility signals a branch-off to a mass-deformed KS black hole in a phase diagram of solutions of 
the cascading theory. However, if = 0, the KT black hole is still unstable but there is no 



associated KS black hole. The chiral limit (3.4) suggests that in this black hole background the 



"chiral tachyons" condensate with finite momentum: the resulting merging phase (and eventual 
endpoint state in a time evolution) cannot be a homogeneous and isotropic geometry described 
by what would be the KS black hole. 

3.2 The mass-deformed Klebanov-Strassler black hole 

In this section we will briefly review the numerical construction of the mass-deformed KS black 
hole (40l . This black hole describes the high-temperature phase of the mass-deformed KS theory. 
The zero temperature phase of this theory is dual to the Kuperstein-Sonnenschein perturbative 



solution 41 when the gaugino masses are equal and to the more general perturbative solution 
of [6 14 when the masses are arbitrary. The infrared expansion of the zero-temperature solution 
is included in the class of perturbative IR solutions constructed in [8,42 . Our purpose is to 
calculate the D3-brane charge at the horizon of the black hole in this theory and to see whether 
the sign of the black hole charge could somehow be made negative by scanning through the 
possible ranges of gaugino masses. 

We are interested in the most general field ansatz of the cascading gauge theory that de- 
scribes its homogeneous and isotropic states, both at zero and non-zero temperature. Such 
ansatz is tailored to search for the (mass-deformed) KS black holes and naturally describes also 
the supersymmetric Klebanov-Strassler (KS) warped deformed conifold solution |1|. It also in- 
cludes as special cases the KT black hole solution and the KT singular supersymmetric solution 
discussed in the previous section, which further obey the constraints imposed by requiring the 
Z2M chiral and U(1)b symmetries. In the Einstein-frame this ansatz is 

ds 2 w = g^(y)dy»dy u + ds 2 Thl , ds 2 Thl = n\{y)gl + Slfty) [gj + gj] + n 2 3 (y) [g\ + g 2 ] , (3.5) 
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for the gravitational field (y denotes the coordinates of M.^ with greek indices //, v = 0, • • • ,4), 
and 

B 2 = hx{y) giAg 2 + h(y) g 3 Ag 4 , 

P3 = 05 A 92, A 04 + ^(y) (0i A 02 - 03 A 4 ) A 05 + (01 A 03 + 02 A 4 ) A d (h 2 {y)) , 



9 



ft 



4Q + h 2 (y) (h 3 (y) - h x {y)) + -Ph^y) 



05 A 03 A 04 A 0i A 02 



(3.6) 



for the fluxes H 3 = dB 2 , F 3 , F 5 , and dilaton $. P is again an integer corresponding to the RR 
3-form flux on the compact 3-cycle (and to the number of fractional branes on the conifold), 



as given in (2.4) 



To make closer contact with the original study of the KS warped deformed conifold [T] we 
use the 1-forms on T 1 ' 1 , 



0i 



a 1 — a 3 



V2 



02 



a 2 — a 4 



V2 



03 



a 1 + a 3 



V2 ' 



04 



a 2 + a 4 



V2 



5 = a 5 , (3.7) 



where 



a 5 = dip + cos Q\d(\>\ + cos 6 2 d<p 2 
a 3 = cos ip sin 9 2 d<fi 2 — sin ipdQ 2 , 



a = — sin9id(pi , a = d9\ 
a 4 = sin ip sin 9 2 d(p 2 + cosipd0 2 . 



(3.8) 



instead of those in (2.3). 



With this ansatz, a Kaluza-Klein reduction of the type IIB action to five dimensions yields 
the effective action |40|, 



S R 



108 



167TG, 



/ vou n^pl (i?i - 1 (v$r 



(V/i 2 ) 2 + 



ttjtt 



h 2 - 



P 



■>()■ \ '" 2 g 



+ 



-2^m{ ia, ' +h2<h3 ~ hi)+1 9 Phi 

where Rio is given by in terms of the 5-dimensional Ricci scalar of the metric g^ u as 



(3.9) 



R 



10 



R 5 + 



+ 



n 2 



2D 5 ln (Oi^^O^) 



( V In fii) 2 + 2 (V In tt 2 ) 2 + 2 (V In fi 3 ) 2 + ( V In (ftifi 2 Q 2 ) ) 



(3.10) 



and V and Ds denote the covariant derivative and the Laplacian in A^s, respectively. This 



action (3.9) also describes the theory governed by (2.5) as a special limit 40 



The general 5-dimensional background geometry with homogeneous and isotropic (but not 



necessary Lorentz-invariant) asymptotic boundary takes the form written in (2.8). It is tailored 
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to search for the (mass-deformed) KS black hole with horizon located at x = 1 and asymptotic 
KT background at x — > 0. It is also useful to encode all the information convening the warp- 
squashed T 1,1 factors ^2,3 (x), the flux-form functions hi j2 ,3 (x) and the dilaton in the 
following set of new functions {Ki, K 2 , K 3 , f a , ft,, f c , h, g} defined as: 



CI; 



H(x) = 

Ui(x) = 



[l = 

(2x 



1,2,3), 



x 



h(x) 



X 



3(2x - x 2 ) 1 / 4 



X 



h\{x) 



1 

p 



12 



Kx{x) -360 ( 



V6(2x - x 2 )V4 

p 



h 3 {x) 



uj 3 (x) 

1 

p 



y/Q(2x - x 2 ) 1 / 4 



12 



K 3 (x) 



36fi 
(3.11 



From the equations of motion of (3.9) we find that G xx (x) can be expressed by an algebraic 
relation that is a function of {K±, K2, K 3 , f a , f c , h, g} and their first derivatives. The non- 
trivial equations of motion can then be written as a system of eight second order non-linear 
coupled ODEs for {Ki, K 2 , K 3 , f a , f b , f c , h, g}. Each of these eight equations is second order. 
It follows that the total differential order of the system is 16. The constant Qq, introduced in 
(3.6) does not appear in the equations of motion. 

The (mass-deformed) KS black hole must solve these equations of motion subject to appro- 
priate boundary conditions at the horizon, x — > 1, and at the asymptotic boundary x — > 0. 

Near the horizon (x — > 1), the equations of motion have the IR series expansion: 



x 



,2n 



1,2,3, 



n=0 

00 

f a = a J2fL (l~*) 2n 

n=0 

00 00 

J2h h n (i-x) 2n , g = goYl 

n=0 



a 



a,b,c, 



(3.12) 



h 



9 k n 



X) 



\2n 



n=0 



and regularity of the expansion in the ingoing Eddington-Finkelstein coordinates requires the 
leading terms of this expansion to be just regular. There are 9 independent IR parameters that 
we choose to be: 



IR independent parameters (9): {k lh0 , k 2h o , faho , faho , fahi , fbha , fcho , h h0 , g h o} ■ (3.13) 

At the UV asymptotic boundary we want the (mass-deformed) KS black hole to asymptote 
to the KT solution [2], while allowing also for the mass-deformation parameters. The latter is 
the leading order contribution of a power series expansion in x and In (a;) of the equations of 
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motion: 



11 °° 
K x =4h a 2 - -P 2 g - -P 2 g Inx + J^Yl k ^ k ^ ln * 
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n=l k 
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n/4 ln k x. 
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^o-^ lnx + ^T]T /i nfc x^ 4 ln fc x, 5 = ^ ( 1 + ^ E 



S'nfc 



n/2 j n fc 



,r 



n=2 fc 



n=2 fc 



(3.14) 



There are 6 independent microscopic parameters. Four of them describe the KT solution, 
namely {P, g , a , h }, and have exactly the physical interpretation already described below 
(3.14). In particular, h can be traded for the dimensionless parameter k s , defined in (2.13), 
which defines the dynamical scale A of the cascading theory. Moreover, the straightforward ex- 
tension (to accommodate for the field content of the KS ansatz) of the three scaling symmetries 
discussed in (2. 14)- (2. 15) still allows to set 

P = g = a = 1 . 



(3.15) 



The other 2 independent microscopic parameters are {kno, faio} and are related to the cou- 
plings of the two dimension-3 operators that explicitly break the chiral symmetry. The explicit 
relation between {k uo , f a io} and these two mass-deformation parameters {/xi, /i 2 } of the cas- 
cading gauge theory introduced in (3.2) is 

faio = (Mi + 4/i 2 k a ) e- k °' 2 , (k no + l -P 2 g a fai^J (SP 2 g + Shoal)' 1 = ^2 e~ ks/2 . (3.16) 

In addition to these microscopic UV parameters, there are 7 extra independent parameters 
{/a30 , ^230 , faio , 9io , fa60 , faio , faso} associated to the VEV's of the operators dual to the fields 
we are solving for. Summarizing, after using the scaling symmetries of our system, we are left 
with 10 UV independent parameters, 

UV independent parameters (10): { k s , k no , f aW , f a30 , k 230 , f a w , 9io , faeo , faio , faso}- 

(3-17) 

The total number of IR and UV independent parameters is 9 + 10 = 19 and is given by (3.13 ) 
and (3.17). On the other hand the equations of motion are a system of total differential order 
16. Therefore, a general mass-deformed KS black hole is a 3-parameter family of solutions. 
We naturally take these 3 parameters to be the microscopic parameters of the cascading gauge 
theory, namely {k s , kno, faio}- We can then use (2.13) and (3.16) to express the results in 
terms of the temperature T/A and the mass-deformation parameters {fix, (12}- Of course, if we 
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Figure 5: Mobile D3-brane charge Q^ 3 for mass-deformed KS black holes with mass- 
deformation parameters f aW = 0.1 and fc 110 = 0, as a function of k s (related to T/A). 



are just interested in the KS black hole that would describe the spontaneously symmetry broken 
deconfmed phase, we set \i\ = /i 2 = and we have simply a 1-parameter family of solutions. 
As emphasized previously, the numerical code does find generic mass-deformed KS black holes 
(whose properties will be discussed below), but when we crank down the mass-deformations to 
zero, no solution with 0\ ^ is found. 

The numerical construction of the mass-deformed KS black holes is done using a shoot- 
ing method similar to the one used in the construction of the KT black hole. In this shooting 
method, the IR power series solution (3.12 ) is constructed up to order n = 1 (inclusive), and the 
asymptotic UV series solution (3.14) is constructed up to eighth order {n = 8, inclusive). One 



then evolves these solutions away from these critical points and do the matching in the overlap- 
ping intermediate region (at x = 0.5). Fixing the microscopic parameters {k s , kuo, f a io}, we 
still have 16 independent IR/UV parameters that are determined by the 16 conditions that fol- 
low from requiring that the 8 functions {K 1 , K 2 , K 3 , f a , f b , f c , h, g} and their first derivatives 
are continuous at x = 0.5. The process is repeated for different values of {k s , k 110 , f a io}- 



3.3 Properties of the mass-deformed Klebanov-Strassler black hole 

Two quantities of utmost interest for our discussion are the D5-brane and D3-brane charges of 



the mass-deformed KS black hole. The former is simply given by (2.22), Qfj£ x = M. On the 



other hand, we can compute the Maxwell D3 charge of the black hole (2.25) using (3.12) to get 



OF 



ft 10 



\ z ^20) i ^20 ^30 



547T 



(3.18) 



Note that the constant constant flo introduced in (3.6) does not appear in the equations of 
motion neither in the flux F§ but it does appear in the contribution B 2 A F% to the D3-Page 



charge (2.23). Below we present the properties of the mass-deformed KS black hole solution. 
As described above this is a 3-parameter family of solutions parametrized by {k s , k lw , f a io}, 
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Figure 6: Left Panel: Mobile D3-brane charge 3 for mass-deformed KS black holes 
with parameters k s = and / a io = {— 0.4, • • • ,0.4}, as a function of the mass-deformation 
parameter k lw . The central red curve/dots with minimum at k 110 = is for f al0 = 0. As 
faio grows increasingly positive, the minimum shifts to the left (black curves are for f aW = 
{0.1, ••• ,0.4}). As f a iQ grows increasingly negative, the minimum shifts to the right (blue 
curves are for f a io = {—0.4, • • • , —0.1}). Right Panel: Minima of the black hole D3-brane 
charge Q® 3 for mass-deformed KS black holes with respect to kuo for k s = 0, as a function of 
the mass-deformation parameter f aW . 
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Figure 7: Same as Fig. [6j but with k s = —0.5. 

the values of which we give as an input to the code. The values of {k^ , k^Q, k^} are then 
determined numerically. For the several families of values of {k s , ki W , f a w} that we have tried, 
it is always true that Q® 3 > 0. 

We outline now our searches. In Fig. [5] we present values of Q® 3 with mass-deformation 
parameters f a ±o = 0.1 and fcno = 0, as a function of k s (related to T/A). Note that smaller 
values of k s bring down the Maxwell D3-brane charge at the horizon — it is thus desirable to 
study as small values of k s as possible. The technical reason that restricts this is rooted in 
the instability of the finite-difference code^jused to simulate mass-deformed KS black holes for 
small k a . The insert in Fig. [5] shows the smallest values of k s we can trust with our current 
numerics. A quadratic fit (in k s ) to these 8 points shows that 



mm 



Qb 3 [ faio = 0.1, fcno = 0; K 



0.000323195. (3.19) 



8 It is possible that these instabilities can be alleviated with spectral methods. We hope to report on this in 
future work. 
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Figure 8: Relation between the mass-deformation parameters f a io and fcno for k s = (Left 
Panel) and k s = —0.5 (Right Panel) at minQ^ 3 . The solid green lines represent the best fit 
linear relation f aW oc ku . 



Figs. [6] and [7] present the results of a different search strategy. Namely, we fix k s (set k s = 
for Fig. [6] and k s = —0.5 for Fig. [7]), and generate families of mass-deformed KS black holes by 
scanning the mass-deformation parameter k uo for discrete set of f aW = {—0.4, ■ • • 0.4} (see the 
Left panels). For each of these curves we use quadratic extrapolation to determine min kll0 [Q® 3 ] 
(see the Right panels). Finally, Fig. M indicate functional relations between non-normalizable 
modes f a io and k U o at minima of for k s = (see the Left panel) and k s = —0.5 (see the 
Right panel). The solid green lines represent the best linear fit kno oc f a io- Note that (3.16) 
implies 

kuo = faio A ^~\ • (3-20) 

Thus, a linear relation between fc 110 and / a io, if true, would imply a linear relation between mi 
and m 2 . From the fits we determine: 



m 2 
mi 



0.072 , for k, = ; 



0.064 , for fc, 



-0.5. 



(3.21) 



Assuming that discrepancy between ^ for different k s is a numerical artifact, an interesting 
question remains as to physical significance of this ratio. 

Although we cannot exclude the possibility that there is some region in the parameter space 
that yields Q® 3 < 0, we could not find it. We take this, together with the fact that the KT 
black hole of the previous section always has positive Q® 3 (recall that for this solution we could 
run the full parameter space since it is a one-parameter solution), as strong evidence for the 
non-existence of KS black holes with Q^ 3 < 0. To summarize, the mass-deformed KS black 
holes we have examined always have positive Maxwell D3-brane charge at the horizon; there 
appear to be no such black holes that have a Maxwell charge whose sign is opposite to that of 
the asymptotic charge of the solution. 



4 Conclusion 

The idea of a Landscape of deSitter vacua in String Theory is founded on two pillars. The 
first is the construction of a very large number of string theory compactifications with stabi- 
lized moduli, either via the KKLT construction [9] or by other mechanism 12 , which always 
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give anti-deSitter vacua. The second is the uplift of these vacua to deSitter, and the most 
common way to do this is to trap anti-D3-branes in warped-deformed conifold-like regions of 
the compactification manifold |9|. The recent analysis of |3]^6j[8] shows that the geometries 
corresponding to these antibranes are singular, and that furthermore this singularity does not 



appeared to be resolved a la Polchinski-Strassler 22 , which is its most obvious resolution 
channel [?] ; this supports the idea that this singularity is not physical, and should therefore be 
discarded. 

In this paper we have looked at the issue from a different perspective: If the singularities 
coming from anti-D3-branes at the bottom of the Klebanov-Strassler throat are physical, then, 



following 28 , it should be possible to cloak them with a regular Schwarzschild horizon. The 
resulting nonextremal geometries would be Klebanov-Tseytlin or Klebanov-Strassler black holes 
that have negative D3 brane Maxwell charge Qf 3 at the horizon. These black holes were 



constructed numerically in 36,37, 40 and we did an extensive scan of all these numerical 
solutions but were unable to find any solution with Q^ 3 < 0. While not a rigorous proof, our 
negative result strongly supports the idea that singularities generated by anti-D3-branes at the 
bottom of the warped deformed conifold are unphysical. This in turn suggests that the most 
common mechanism for uplifting the landscape of AdS vacua one obtains from string theory 
compactifications with stabilized moduli to deSitter does not work, and hence string theory 
may only have a landscape of AdS vacua but not a landscape of deSitter vacua. 

Furthermore, we also found that the charge of a given KT or mass-deformed-KS black hole 
is not an independent parameter, but is completely determined by the temperature and the 
gaugino masses. This implies that if one is to perform a gedanken experiment that consists 
of lowering an anti-D3 brane into this black hole keeping the temperature fixed, the charge of 
this black hole will become a bit smaller for a moment, but then will immediately go back to 
its previous value by absorbing charge from the surrounding fluxes. Hence, this black hole acts 
as a catalyst for brane-charge annihilation for arbitrarily small values of the temperature, and 
this may suggest that adding an antibrane to a vacuum Klebanov-Strassler solution will cause 
the surrounding flux to immediately annihilate it [45]. 

It would be interesting to further push the limits of the KT/KS black hole parameter 
space (driven by the development of more efficient numerical techniques) in order to establish 
more completely that regular nonextremal geometries with Q^ 3 < do not exist. It would 
also be interesting to compute D3 brane charge of the bottom of de-Sitter deformed KT/KS 



geometries 43 44 . We hope to report on this in future work. 
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